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ABSTRACT - This report provides information 
on the current status of NBS work on mea- 
surement technology for semiconductor mate- 
rials^ process control, and devices. Empha- 
sis is placed on silicon and silicon-based 
devices. Highlighted activities include 
semi -insulating C^As characterization, an 
automatic scanning spectroscopic ellipsome- 
ter, linewidth measurement and coherence, 
bandgap narrovdng effects in silicon, the 
evaluation of electrical linewidth uniform- 
ity, and arsenic-implanted profiles in sili- 
con, In addition, recent publications and 
publications in press are listed. The re- 
port is not meant to lie exhaustive; con- 
tacts for obtaining further information are 
listed, 

KEY WORDS - Electronics; Integrated cir- 
cuits; measurement technology; microelec- 
tronics; semiconductor devices; semiconduc- 
tor materials; semiconductor process con- 
trol; silicon. 


Preface 


This report covers results of work during the fifty-second quarter of the NBS Semiconductor 
Technologi' Program, This Program serves to focus NBS research on improved measurement 
technology for the use of the semiconductor device community in specifying materials, 
equipment, and devices in national and international commerce, and in monitoring and con- 
trolling device fabrication and assembly. This research leads to carefully evaluated, 
well-documented test procedures and associated technology which, when applied by the indus- 
try, are expected to contribute to higher yields, lower cost, and higher reliability of 
semiconductor devices and to provide a basis for controlled improvements in fabrication 
processes and device performance. By providing a common basis for the purchase specifica- 
tions of government agencies, improvsd measurement technology also leads to greater economy 
in government procurement. Financial support of the Program is provided by a variety of 
Federal agencies. The sponsor of each technical project is identified at the end of each 
entry in accordance with the following code: 1, The National Bureau of Standards; 2. The 
Defense Advanced Research Projects Agency; 3, The Division of ElectrLc Energy Systems, 
Department of Energy; 5, The Defense Nuclear Agency; 6, The C, S, Draper Laboratory; 10, 
The Naval Weapons Support Center; 11, The Solar Energy Research Institute; 13, The Lewis 
Research Center, National Aeronautics and Space Administration; and 15, The Naval Ocean 
Systems Center, 

This report is provided to disseminate results rapidly to the semiconductor community. It 
is not meant to be complete; in particular, references to prior work either at NBS or else- 
where are omitted. The Program is a continuing one; the results and conclusions reported 
herein are subject to modification and refinement. Further information may be obtained by 
referring to more formal technical publications or directly from responsible staff members, 
telephone: (301) 921-listed extension. General information, past issues of progress 
briefs, and a list of publications may be obtained from the Semiconductor Materials and 
Processes Division, National Bureau of Standards, Washington, D.C, 20234, telephone: 
(301) 921-3786. 
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Bandgap Narrowing Effects in Silicon 

The effects of heavy doping on the den- 
sity of states and the energy bands are 
important for the interpretation of 
device behavior. The average shift of 
the conduction and valence bands due to 
the interaction of electrons and holes 
with donor impurities has been calcu- 
lated by using a partial wave analysis. 
In order to demonstrate feasibility, or 
proof of concept, this work was re- 
stricted initially to states near the 
band edge since only the i = 0 partial 
wave would be important. The model 
Yukawa potential was also approximated 
with a potential due to Bargmann for 
which the S, = 0 phase shift for the 
latter potential can, unlike that for 
the Yukawa potential, be expressed as a 
function of the wavenumber in closed 
form. This decreased the computational 
expense considerably. Tlie success of 
the above method justified the following 
improvements. The range of validity in 
wavenumber was extended by including the 
partial waves with Jl = 1 and = 2 as 
well as A = 0. The Yukawa potential was 
also used throughout since numerical 
solutions for the phase shifts are nec- 
essary whenever £ > 1 even for the po- 
tential due to Bargmann. An additional 
refinement has been the self-consistent 
calculation of the screening length 
associated with the Yukawa potential by 
satisfying the Priedel sum rule. The 
results of this calculation show impor- 
tant differences between this technique 
and the Fermi -Thomas approach. 

The shifts in the energy gap resulting 
from these calculations are smaller than 
the shifts measured optically. The 
difference between the experimental data 
and these results is due to many body 
effects, primarily exchange and correla- 
tion. In fact, the results of Mahan for 
these contributions, which are valid at 


0 K, when added to these results yield 
good agreement with the data at 0 K. At 
room temperature, it is not possible to 
make such a direct comparison. Good 
agreement also is obtained between these 
results for the Fermi energy measured 
with respect to the conduction band edge 
and the corresponding optical data. The 
bandgaps from optical data are much 
smaller than the bandgaps from electron 
device experiments. The implication is 
that other effects, such as variation of 
the recombination lifetime, should be 
included in the analysis of the electri- 
cal measurements. 

An important result of these calcula- 
tions is that not only is the energy gap 
narrowed, but the shape of the bands is 
changed as well. If one were to try to 
fit parabolic bands to the perturbed 
bands over a given region necessary for 
co ^iutations, the effective mass of the 
new parabolic bands would be different 
from the unperturbed values and would 
depend upon the extent of the region 
fitted. 

These results have been applied to pro- 
cess modeling by considering the effect 
of bandgap narrowing on the diffusion 
coefficients of donor ions at process- 
ing temperatures. These coefficients 
depend strongly on the ratio of the 
electron density to the intrinsic carri- 
er density. An important question is 
whether bandgap narrowing effects sig- 
nificantly alter the interpretation of 
diffusion data as would be implied by 
the results of other workers for donor 
and vacancy activity coefficients. The 
result of applying this model to a 
sample containing a donor density of 
1.0 X 10^® cm”^ at 1100°C is a bandgap 
reduction, -AEg, of 123 meV. However, 
the ratio of the effective intrinsic 
carrier density, n^e to the undoped 
value, n^, is only 1.15, which is much 
less than one obtains from the expres- 
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sion exp(-ABg/2kT), l*e., 1 •68> The 

reason for this large difference is that 
the density of states is highly dis- 
torted from its parabolic form« The 
model described here is very different 
from previous models and in general 
yields much smaller values of the ratio 
nie/ni. [Sponsor: 11 

(J. R. Iiowney and H. S. Bennett; x3625) 

Arsenic-Implanted Profiles in Siicon 

Device performance will become more 
sensitive to the three-dimensional na- 
ture of implanted and annealed distribu- 
tion of the dopants with the trend to- 
ward smaller device dimensions for very 
high speed, very large scale integrated 
circuitst This may serve as the limit- 
ing factor in device performance. The 
two-dimensional distribution of arsenic 
implanted into silicon was studied in 
this work. 

The basis of the present understanding 
of the vertical distribution of ions 
implanted into an amorphous target is 
the Winterbon code which is based upon 
he numerical solution of the Boltzmann 
equation for the moments of the distri- 
bution. Extension of these calculations 
to two dimensions is nontrivial. How- 
ever, there are Monte Carlo codes, spe- 
cifically the TRIM (Transport of Ions 
into Matter) due to Biersack and Hagg- 
mark, which may be used to construct the 
two-dimensional distribution of ions 
implanted into amorphous targets. Both 
of these codes have been used to con- 
struct the vertical distribution of 
arsenic implanted into silicon at 150 
keV and yield similar results for the 
same choice of the ion-atom potential. 
In addition, these results have been 
compared with experimental data from 
SIMS and the three are in good agreement 
with one another. Further, the TRIM 
code has been used to construct the two- 
dimensional distribution of arsenic 
implants and these results serve ast the 
basis for a physical understanding of 
the experimental results discussed 
here. 


In addition, the predictions of the 
SUPREM-II code for process modeling 
supplemented with the more recent under- 
standing of the arsenic diffusion pro- 
cess have been employed in the evalua- 
tion of the results of annealing studies 
of arsenic- implanted regions. 

In the experimental study, specially 
prepared specimens fabricated at Hughes 
Research Laboratories ( were employed. 
Specifically, boron-doped silicon wafers 
were prepared with a grating pattern 
structure on the surface. This pattern 
served to mask out the Incident arsenic 
ions during the implantation process 
allowing the arsenic to enter only the 
silicon lattice in parallel stripe re- 
gions. The wafers were annealed under 
various time -temperature conditions in 
order to ascertain the effects of the 
anneal cycle on the distribution of the 
arsenic atoms due to both the original 
implantation process and the anneal- 
induced redistribution. 

After the patterned wafers had been 
implanted and annealed, the photoresist 
and aluminum stripe pattern were removed 
and the wafers wore cleaned with stan- 
dard peroxide solutions. Finally, the 
wafers were cleaved perpendicular to the 
grating fingers. Two anneal tempera- 
tures were used in the study: 600**C, to 

remove the initial implant damage and to 
activate the arsenic with minimal redis- 
tribution? and 1000®C, to determine the 
effect of the anneal-induced diffusion 
on the final junction location (both 
vertical and lateral). To determine the 
junction location, both etch and EBIC 
(Electron Beam Induced Current) tech- 
niques were used. It is important to 
note that this is the first time that 
these techniques have been used in the 
depth scale required for VLSI and VHSIC 
structure investigation. 

The etch technique, even though there 
were some problems with the application 
of "wet chemistry" to submicrometer 
dimensions, showed that the 600®C an- 
nealed arsenic had no appreciable signs 
of lateral channeling, while the 1000“C 
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annealed arsenic had significant signs 
of lateral diffusion. 

The EBic data are in substantial agree- 
ment with this conclusion. In princi- 
ple, the EBIC data should be able to 
provide not only the location of the 
vertical and lateral junction but also 
the distribution of activated dopants. 
However, this requires the numerical 
solution of the semiconductor transport 
equations in the presence of electron- 
hole pair generation due to the electron 
beam. 

The above observations are in substan- 
tive agreement with the results of the 
two-dimensional Monte Carlo (TRIM) cal- 
culations in that they show very little 
lateral spreading of arsenic when im- 
planted into silicon. Ihis can be un- 
derstood on the grounds that arsenic is 
more massive than silicon and is unlike- 
ly to be deflected in the lateral direc- 
tion by the silicon atoms. On the other 
hand, dopants lighter than silicon, 
boron, for example, might show appreci- 
able lateral spreading. the diffusion 
predictions of the augmented SUPREM-II 
code agrees well with the vertical junc- 
tion location. [Sponsor: 15] 

(J. Albers, J. R. Ehrstein, 
D. R. Myers,* and P. Roitman, x3625) 

Electrical Linewidth Uniformity Evahiation 

An electrical measurement method using 
the cross-bridge test structure to eval- 
uate linewidth variation associated with 
integrated circuit lithography has been 
developed. Arrays of identical cross- 
bridge resistors were patterned in an 
aluminum layer on a silicon wafer using 
a single photomask. Electrical test 
results from these structures were used 
to measure the uniformity of the line- 
width across a wafer and across a chip. 
The data were analyzed to evaluate the 
linewidth uniformity of a step-and- 
repeat system. Techniques were devel- 
oped to identify and remove systematic 
linewidth variations introduced by other 

• Present: address: Sandia Laboratories, Albuquerque, 

New Mexico# 


aspects of the photolithography process 
and to determine the precision of the 
measurement. 

Analysis was made from data of test 
chip MBS- 23. NBS-23 is a test chip 

consisting of an 20 by 12 array of iden- 
tical cross bridges with design line- 
widths of 6 ym. Hie array fits within a 
square 250 mil (6.35 mm) on a side and 
is repeated across the photomask. A 
wafer is then fabricated by delineating 
the cross bridges (using contact print- 
ing and wet chemical etching) in an dOO- 
nm thick aluminum film evaporated on an 
oxidized silicon wafer. Hie automated 
dc test system used to make the measure- 
ments consists of a minicomputer- 
controlled wafer prober, scanners, digi- 
tal current source, and digital volt- 
meter. Approximately 3 s were required 
to measure each cross bridge. 

In order to determine linewidth uniform- 
ity with this test chip, it is Important 
to know the measurement precision which 
determines the achievable resolution. 
Linewidth measurement precision ^is lim- 
ited by the test system and the test 
structure. Hie precision of the teat 
system was determined by repeated mea- 
surements of the electrical linewidth on 
the same cross bridge and a computation 
of the sample standard deviation. The 
precision was found to be 0. 002 nm ( 1® ) 
for sixty measurements. To determine 
the usable precision of the test struc- 
ture, a sample of the measurement data 
was analyzed. Test results from one row 
of cross bridges across the wafer repre- 
senting seven adjacent test chips, each 
chip containing 20 cross bridges, are 
shown in the accompanying figure. Three 
types of variations can be seen: a 

random variation, a systematic intrachip 
variation, and a systematic Intrawafer 
variation. The magnitude of the random 
variation determines the achievable 
measurement precision of the test struc- 
ture. To estimate the magnitude of the 
random variations correctly, the system- 
atic variations in the data must be 
removed. This was accomplished by fit- 
ting the data to a fourth-degree poly- 
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NBS-23 data of linewidth measurements 
from one row across the center of the 
wafer (0) and the fitted function. 


and 8Ubstrate» respectively^ and is 
the optical phase difference at the line 
edge which occurs on reflectance (or 
transmittance). It has now been shown 
that this equation is a special case of 
a more general expression for determin- 
ing the threshold corresponding to the 
line edge location which can be applied 
to any state of partial coherence and 
allows the inclusion of known optical 
aberirations such as defocus. In the 
more general case, the threshold Tpc 
is defined as the ratio of the optical 
intensity in the image at the line edge 
to the reflected (or transmitted) opti- 
cal intensity from the substrate, i.e., 
using a Fourier series representation of 
the image, 


nomial which estimated the intrawafer 
and intrachip variation. The result of 
this fit can also be seen in the figure 
in the form of the connected line vrhlch 
is the sum of the two estimated mathe- 
matical functions. The residual stan- 
dard deviation after subtracting the 
estimated functions from the original 
data is 0.03 Mm (la). 

Electrical linewidth measurements with 
the cross-bridge test structures fabri- 
cated on silicon wafers can resolve 
variations in linewidth to a precision 
of 0.03 Mm (1o). ITiis level of preci- 
sion allows this technique to be used as 
a valuable tool for evaluating the per- 
formance of lithography systems capable 
of submicrometer resolution. [Sponsor: 
1] (D. Yen and L. W. Linholm, x3541) 
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coherent Illumination, the Fourier coef- 
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Linewidth Measurement and Coherence 

Earlier work on optical methods for 
linewidth measurement for IC photomasks 
and wafers led to the development of the 
method of coherent edge detection. In 
the coherent case, the optical intensity 
or threshold corresponding to edge loca- 
tion was given by 



where Iq and 1^ are the reflectances 
(or transinittances) of the etched layer 


+ A. 


n-n 



t 


where 



*o "5 I' ' 

This equation assumes well-resolved 
lines where, for convenience only, the 
linewidth has been taken to be half the 
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period (W/P = 1/2)* The funotion T is 
called the transmission cross coeffi- 
cient and characterizes the 2-D (two- 
dimensional) optical system includino, 
the partial coherence in the illumina- 
tion. For a 1-D (one-dimensional) line 
object 

00 00 

h) ” f f A(eMi")K?i + h") 


X f*(?2 + h") d?" dn" , 

where A (?"/ h") is the 2-D intensity 

distribution in the condenser aperture 
(wide-field Kohler or critical illumina- 
tion assumed). The funotion A is some- 
times called the effective source (as- 
signed to be incoherent; i.e., the coher- 
ence interval is small compared to the 
effective source size). P is the pupil 
function of the optical system in nor- 
malized coordinates. 

Even with spatially incoherent sources, 
a finite degree of partial coherence is 
generated with propagation of light 
through the optical system due to the 
finite dimensions of the source and 
apertures. For wide-field microscope 
illumination, and either critical or 
Kohler illumination, the degree of co- 
herence (spatial) is determined by the 
ratio of condenser to objective numeri- 
cal aperture, frequently called the 
coherence parameter R. The transmission 
cross coefficients represent a correla- 
tion of these aperture functions. In 
the first figure, Tpc is plotted for 
matched condenser and objective numeri- 
cal apertures. Imaging with matched 
apertures was once thought to be "inco- 
herent." However, the obvious sensitiv- 
ity of the edge detection threshold to 
the optical phase difference at the line 
edge is further proof that this is not 
the case, and the matched aperture con- 
dition instead represents an intermedi- 
ate case of partially coherent imaging. 
The 33% threshold used for edge detec- 
tion with matched apertures was first 
noted by Charman. As may be seen in the 
figure, a 33% threshold applies only to 



Threshold TpQ corresponding to edge 
location for coherence parameter R = 2 
as a function of Tq (IQ/Iff), 
Curve-to~curve parameter is optical 
phase difference a) it, b) n/2, c) 

0, The sot incoherent edge detection 
threshold is shown (curve d) for compar- 
ison, 

high contrast ® 

jects where the optical phase difference 

may be Ignored. 

Although it is not possible to achieve 
effectively incoherent illumination in 
high resolution wide-field microscopy 
(with a dry objective, 0.9 N.A., R val- 
ues appreciably greater than 1.0 are not 
possible), it is of interest to many 
microscopists to determine the value of 
R required for effectively incoherent 
illumination. The second figure (next 
page) shows the threshold variation for 
R = 2. This figure Indicates that a 
phase difference of n, which should not 
be visible in incoherent illumination, 
for R = 2 would image as a dark line 
with a modulation of 12% and thus would 
be visible to the eye as well as .to a 
photodetector. Hierefore, effective in- 
coherence for sharp edges requires an R 
value greater than 2. 

The formalism of the transmission cross- 
coefficients allows the inclusion of 
optical aberrations such as defocus in 
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Threshold TpQ cor res pondin';^ to edge 
location for coherence paxaeeter R “ 2 
as a function of Tq 
C urve-to-curve parameter is optical 
phase difference 4’c* 

0, The 50i incoherent edge detection 
threshold is shown (curve d) for compar- 
ison. 



Variation of TpQ with defocas param- 
eter a^ for Tq w 0.1, for different R 
values t dashed lines. R - 1/5, solid 

lines R « 2/3, The coherent edge detec- 
tion threshold Tq is indicated in each 
case and is equivalent to the value of 
fpc zero defocus and R - 1/5. 
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terms of the pupil funation of the imag- 
ing optics. For defocus/ the pupil 
function liecomes 

%focu|P> " P - (5^ + n^) 

Where a^. indicates the amount of defoous 
in waves (units of wavelength) and )« is 
the wave number. The third figure shows 
the variation of the edge detection 
threshold with focus for selected values 
of Tq and Note that with in- 

creased coherence the edge detection 
threshold is less sensitive to defooUs 
and thus measurement precision improves. 
[Sponsors 1) (D. Nyyssonen, X3621 ) 

Semi-Insulating GaAs Characterization 

A measurement technique for charac- 
terizing electrically active traps in 
semi -insulating (Si) semiconductors has 
shown potential for further development. 
The technique, termed photoinduced 
resistance-deep level transient spec- 
troscopy (PR-DLTS), permits electrical 
characterization of high resistivity 
semiconductors such as Rl GaAs, for 
which many conventional electrical tech- 
niques are inadequate. In PR-DLTS, 
monochromatic light is used to create 
excess free carriers in the material 
which are then trapped by the deep lev- 
els. After the light is extinguished, 
the excess free carriers recombine, and 
there occurs a transient change in the 
specimen conductivity as the carriers, 
trapped at the deep levels, are emitted 
and the system returns to thermodynamic 
equilibrium. The rate at which the 
trapped carriers are emitted ie a char- 
acteristic function of the particular 
trap involved and the temperature. This 
emission rate is measured as a function 
of temperature resulting in a signature 
of the particular trapping center. 

The transient conductance is measured in 
an impedance bridge at 10 MHz and signal 
processing is performed as in conven- 
tional DLTS. In the PR-DLTS measure- 
ment, however, specimen capacitance is 
fixed by a parallel plate test structure 
with the SI GaAs assuming the role of a 
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Photoinduced resistance-DLTS measurement 
of semi- insulating QaAs for LEC-grown 
material (a) and HB-grown material (b). 


dielectric. Light illuminates the SI 
GaAs surface through spaces in an alumi- 
num grid contact. Hy phase- sensitive rf 
detection techniques the resistive por- 
tion of the transient is monitored as 
opposed to the capacitive portion in 
conventional DLTS. Results to date show 
that the method is capable of discrimi- 
nating between variously grown specimens 
of SI GaAs as displayed in the accom- 
panying figure. PR-DLi’S spectra are 
presented for two specimens of SI GaAs 
grown by the same manufacturer, but by 
two different techniques: liquid encap- 
sulated Czochralski (LEO) and horizontal 
Bridgman (HB). As can be seen, four 
peaks are evident in both the LEG and HB 
spectra, each peak corresponding to a 
specific deep level. Interpretation of 
these data in terms of particular impu- 
rities, densities, and energy levels has 
not yet been established. 

PR-DLTS shows promise as a tool to elec- 
trically characterize device starting 
material as well as to correlate materi- 
al properties with device performance 
parameters. It offers the advantage. 


ovor other oloctrioal and photo- 
olQotronic methods, of noncritioal con- 
tacts. Ihe system is also readily con- 
structed by simple modifications of a 
conventional DLTS system* 

Recent findings have shown the PR-DLTS 
spectra to be a function of light exci- 
tation wavelength, intensity, and dura- 
tion. ohese effects are presently being 
examined. In addition, experiments are 
planned to correlate the PR-DLTS data 
with other electrical, optical, and 
physical measurement techniques as well 
as to interpret the results in terms of 
specific deep level states. (Sponsor: 
1] (A. C. Soabaugh, x3625) 

Automatic Spectroscopic Ellipsometer 

A new scanning automatic spectroscopic 
ellipsometer is being built in the Semi- 
conductcr Materials and Processes Divi- 
sion. TTiis ellipsometer will be a vari- 
able angle of incidence instrument that 
can be run in either an automated prin- 
cipal angle of incidence mode, rotating 
analyzer mode, or conventional automatic 
null mode. Each method has its own 
advantages depending upon the user's 
needs, but an important feature will be 
the ability to intercompare these meth- 
ods. Itie mechanical design is one based 
on rotating goniometers powered by step- 
ping motors as used, for example, in x- 
ray diffraction instruments. The goni- 
ometers have very high accuracy and 
precision required for this ellipsome- 
ter. The light sources will be a stan- 
dard HeNe laser (632. 8-nm wavelength) 
and also a HeCd (441.6-nm wavelength) 
laser. In addition, a prism-grating 
monochromator (186- to 2650-nm wave- 
length, corresponding to 6.7 to 0.47 eV) 
will be used for spectroscopic studies. 
The optics will include standard ellip- 
sometric polarizers and compensators. 
The automation will be carried out by 
the use of a computer interfaced to 
address the components of the instru- 
ment, which include the stepping motor 
controlled goniometers and a photomulti- 
plier or solid-state light detector. 
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Communion tiona from roadora having an 
intereat In precise ellipsomatrio mQa~ 
guremants are welcomed. [Sponsor! 1) 

(G» A. Candela ana D. Horowitz, x362S) 


Erratum lr» ASTM F 723-81 


A printing error exists in Astm P 723-81 
for conversion between resistivity and 
dopant density. This standard practice 
is based on NB8 data. The correct equa- 
tion to convert resistivity to dopant 
density for boron-doped silicon is: 


1.330 X 

P_ 


10 


16 


1.0824 X 10 


17 


p( 1 + (54.56 


[Sponsors! 1,2] {J. R. Ehrstein, x3625) 
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